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Motivated by recent experiments with confined binary liquid mixtures near 
their continous demixing phase transition we study the critical behavior of a 
system, which belongs to the Ising universality class, for the film geometry with 
one planar wall chemically structured such that there is a laterally alternating 
adsorption preference for the species of the binary liquid mixture. By means 
of Monte Carlo simulations and finite-size scaling analysis we determine the 
critical Casimir force and the corresponding universal scaling function. 
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1. Introduction 



[~^. ■ If a field exhibiting fluctuations with long wavelengths is confined between 

CD I surfaces, long-ranged effective forces between them arise. The order param- 

eter of a second-order phase transition represents such a field. The resulting 
effective force is known as the critical Casimir force. This phenomena, first 
predicted by Fisher and de Gennes 1 is the analogue of the Casimir effect 
in quantum electrodynamics. 2,3 The critical Casimir force is characterized 
by a universal scaling function, which is determined by the bulk and sur- 
face universality classes (UC) of the confined system. It is independent of 
the microscopic details of the system, and it depends only on a few global 
and general properties, such as the number of components of the order pa- 
rameter, the shape of the confinement, and the boundary conditions (b.c.) 
there. 

In recent years, the critical Casimir effect has attracted numerous ex- 
perimental and theoretical investigations; see Refs. 3,4 for recent reviews. 
Wetting layers of binary liquid mixtures 5 have provided indirect evidence 
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and measurements of critical Casimir forces. Recently, a direct determi- 
nation of the critical Casimir force has been reported 6 using individual 
colloidal particles immersed in a binary liquid mixture close to its critical 
demixing point and exposed to a planar wall. 

Quantitatively reliable calculations of critical Casimir forces for later- 
ally homogeneous b.c. have been obtained recently by means of Monte 
Carlo (MC) simulations. In this context, the Ising UC, 8 which describes 
the demixing transition in a binary liquid mixture, has been investigated 
and the critical Casimir force has been determined. 

Experiments with binary liquid mixtures have also been used to study 
critical Casimir forces acting on a colloid in front of a chemically struc- 
tured substrate, 9 leading to a laterally varying adsorption preference. Such 
a system has been investigated theoretically for planar geometries within 
mean- field theory 10 and for a curved geometry within the Derjaguin ap- 
proximation. 11 Here we consider the Casimir force in the film geometry 
L x iii x Lii, Lu ^$> L, such that the free energy density T per volume 
decomposes as 

?(0,L,L\\) = / bulk (/3) + i/ s „ rf (/3) + i/ cx (/?,i,i||), (1) 

where /buik(/3) and / sur f(/3) are the bulk and the surface free energies per 
volume and area, respectively, in the thermodynamic limit L\\,L — > oo, 
and /ex is the excess free energy. The Casimir force Fc per area and per 

p- 1 = k B T is 



p Wex) 



c)L 



(2) 



According to finite-size scaling, 12 close to a continuous phase transition 
Fq attains the following asymptotic scaling behavior: 

F C = -If? (t(L/^) 1/v ,p) , t = (T- To)/T c , (3) 

where Tc is the bulk critical temperature, p = L/Ln is the aspect ratio, and 
£o" is the non-universal amplitude of the correlation length £ = £rj t~ v in the 
disordered phase. The scaling function 9{x,p) is universal, i.e., it depends 
only on the UC of the bulk phase transition and on the surface UC, which is 
determined by the b.c. applied on the two surfaces. In the case of a confined 
binary liquid mixture, the order parameter <j) is the difference between the 
local and the critical bulk concentration of one of the two components. 
The substrates generically prefer one species over the other, so that \<j>\ is 
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Fig. 1. Film geometry with a single chemical step; S is the spin value at the surface. 

enhanced near the surface: this corresponds to the so-called extraordinary 
or normal surface UC 13 which, in a lattice model, can be implemented by 
external fields acting at the surfaces. For laterally homogeneous substrates, 
one has to distinguish between the case of surfaces with the same (++) and 
opposite (H — ) adsorption preferences. For these systems laterally varying 
adsorption preferences can be realized. 9 

Motivated by corresponding experimental results, 9 we have performed 
MC simulations of a representative spin model on a lattice belonging to the 
Ising UC. 14 As illustrated in Fig. 1, the spins on the upper surface are fixed 
to +1. The lower surface is divided into two halves, one with spins fixed to 
—1 and the other with spins fixed to +1; the lattice constant is set to 1. 
This mimics a single chemical step as the simplest lateral heterogeneity. 

2. Simulation method 

In order to calculate the Casimir force, we use the so-called coupling pa- 
rameter approach. 8 Given two systems with Hamiltonians Hq and H\ and 
the same configuration space, the free energy difference AF between the 
two systems can be computed as 

AF= f d\(AH) x , AH = U X -Ho, (4) 

Jo 

where (. . .)\ denotes the thermal average over the ensemble described by the 

crossover Hamiltonian TL\ = (1 — X)Hq + XH\, with A £ [0, 1]. This average 

can be evaluated by standard MC simulations, followed by numerically 

carrying out the integral in Eq. (4). We have used Eq. (4) by considering 

"Hq as the Hamiltonian of the system shown in Fig. 1 and Hi as that of the 

corresponding film but with thickness L — 1 and an additional, decoupled 

two-dimensional (2D) layer. One obtains 

^ = -/buik(/3) + / 2 d(/3) + -^9 (t(L/^f'\p) , (5) 
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Fig. 2. The amplitude © (p) of the critical Casimir force at T c for the system shown in 
Fig. 1 as a function of the aspect ratio p = L/Lu ■ 



where /2D (/?) is the free energy density of the decoupled 2D layer and we 
have set the lattice constant to 1, so that each term is dimensionless. We 
note that Eq. (5) is correct only up to corrections to scaling. 



3. Results 

At the bulk critical temperature the expression in Eq. (5) reduces to 

1 f 



AF 



Pc L a 



(6) 



where (p) = 8 (0,p) is the critical Casimir amplitude. By fitting the MC 
data to this equation and by adding the expected additional scaling correc- 
tions, one can determine the critical Casimir amplitude. 14 

As shown in Fig. 3 the Casimir amplitude varies linearly for small p. We 
find 6(1/6) = 2.048(6), 6(1/8) = 2.126(5), 6(1/10) = 2.183(6), 6(1/12) = 
2.223(7), and 6(p -> 0) = 2.386(4). In the limit p -> the force is expected 
to be the average of the force for laterally homogeneous (++) and (H — ) b.c. 

At criticality 8 (6 ++ + 6^ )/2 = 2.33(4), in agreement with our results. 

The calculation of the full scaling function 6(x,p) requires the subtraction 
of the L— independent terms in Eq. (5) (see Ref. 14 for details). In Fig. 3 we 
report the resulting scaling function for various aspect ratios. We also show 
a comparison with the scaling function obtained as the average of those for 
homogeneous (++) and (H — ) b.c, based on the data of Ref. 8. For p — > 
we find good agreement with our results. 



July 22, 2010 14:53 



WSPC - Proceedings Trim Size: 9in x 6in paper 



6(x,p) 





1 ' 1 ' 1 


i 


i i i i i 


T# 3K * 


T- T 






. p=l/6 




A 
_A 


* *' 






. p=l/8 




" • 


; 1 * t 






. p=l/10 




• 


• . t * T 






» p=l/12 






i 






T p">0 

* mean(IV) 












- 








~ 


». 




- 


_ 


V 




_ 






¥ 


- 






% , 









i 


*»* » 

i f r i . 



4 8 12 

X =t(L/0 1/v 



16 



20 



Fig. 3. The universal scaling function of the critical Casimir force for the system shown 
in Fig. 1 and for various aspect ratios p = L/Lu. We also show the scaling function 
obtained as the mean value of those for laterally homogeneous (H — h) and (H — ) b.c. based 
on the approximant IV presented in Rcf. 8. Error bars are smaller than the symbol size. 
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